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Conductance of zigzag interfaces between graphene sheet and normal metal is investigated in 
the tight-binding approximation. Boundary conditions, valid for a variety of scattering problems, 
are constructed and applied to the normal metal - graphene - normal metal (NGN) junctions. At 
the Dirac point, the conductance is determined solely by the evanescent modes and is inversely 
proportional to the length of the junction. It is also independent on the interface resistance. Away 
from the Dirac point, the propagating modes' contribution dominates. We also observe that even 
in the junctions with high interface resistance, for certain modes, ideal transmission is possible via 
Fabry- Perot like resonances. 



I. INTRODUCTION 

Recent experimental studies^!" uncovered unusual 
properties of graphene (graphite monolayers and bilay- 
ers), strongly contrasting with the common knowledge 
inherited from studies of metals. This difference origi- 
nates from the fact that electrons in graphene monolayers 
obey the Dirac (rather than the Schrodinger) equation. 
Thus, one has an opportunity to study properties of the 
Dirac fermions in a table-top experiment. Predictions 
of relativistic effects including the Klein tunneling^ and 
Zitterbewegung£ have been made. "Ordinary" phenom- 
ena, e.g., quantum Hall effect, weak localization^ or An- 
dreev reflection^, are also strongly modified in graphene 
as compared to normal metals. 

The most easily accessible measurements in graphene 
are those of electrical transport. Theoretically, one way 
of understanding them is to extend the Landauer theory 
to graphene sheets, considering them as a junction be- 
tween two reservoirs. So far, a common point was to de- 
scribe reservoirs as bulk disordered graphene ^ 10 ' 11 . This 
approach considerably facilitates calculations; however, 
its relation to the experimental situation, with contacts 
made of normal metals, requires additional clarification. 
Indeed, the bulk in the available graphene devices is de- 
fect freei^. Thus, the major source of electron scattering 
in graphene are the boundaries, in particular, contacts, 
and even qualitative understanding of electric transport 
can not be achieved without careful consideration of elec- 
tron behavior in the contacts. 

Experimentally, graphene flakes are contacted by tun- 
nel junctions located on top of the flakes. Conceptually, 
we can discriminate between three types of the junctions. 
One situation is when tunneling from the normal reser- 
voir to graphene occurs just at one point — for instance, 
due to the non-uniform thickness of the oxide layer in 
the contact in combination with the exponential depen- 
dence of the tunneling amplitude of this thickness. Ob- 
viously, in this case the voltage between normal reservoir 
and graphene drops at the junction, and the chemical 



potential in graphene can be considered as fixed. For 
this point-like tunneling the resistance of a NGN junc- 
tion depends on the distance between the two tunneling 
points. 

More interesting, and apparently more experimentally 
relevant is the situation when tunneling occurs at many 
points, covering a large area (Fig. QJ. In this case, there 
is no voltage drop on the contact, i.e. it becomes Ohmic. 
In this area under the contact, the wave functions of 
electrons from the normal reservoir and graphene quasi- 
particles hybridize, forming a new substance, a hybrid 
between graphene and normal metal. The nature of this 
substance depend on the exact properties of the contacts, 
however, it is reasonable to assume that generally it will 
be either closer to graphene, or to a normal metal. In the 
former case, effective description of a NGN system as a 
graphene - graphene - graphene (GGG) contact, where 
the voltage is applied to the graphene, is appropriate. In 
the latter case, the system is essentially a planar NGN 
junction, where the voltage is applied to the normal metal 
and drops at the interface between the normal metal and 
graphene. It is also clear that properties of such contacts 
will considerably depend on the relation between lattice 
periods of graphene and normal metal, for instance, on 
whether the lattices are commensurate or not. 

Below, we discuss the latter type of junctions — pla- 
nar NGN junctions, to complement the earlier studies of 
GGG junctions. A first step in this direction has been re- 
cently made by Schomerus-^ who compared resistances of 
NGN contacts with a zigzag interface and GGG contacts. 
He considered the special case of equal overlap integrals 
between all neighboring sites in the tight-binding model. 
He found that if the graphene sheet is biased to the Dirac 
point, so that there are no propagating modes through 
graphene, the difference between NGN and GGG junc- 
tions is only quantitative. Outside this regime, the be- 
havior of the contacts may strongly depend on the type 
of the leads and the interface. Below, we study this de- 
pendence for the simplest case of square-lattice leads and 
zigzag interfaces, and find a number of interesting results. 

In this Article, we consider within the tight-binding ap- 
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FIG. 1: (a) Typical experimental contact configuration: 
Metallic leads overlaying graphene. The contact occurs over 
many microlinks which leads to equilibration of the chemical 
potentials in the metal and in graphene. In the case of tight 
coupling under the contact the band structures of graphene 
and metal are expected to "fuse", (b) Theoretical model of 
a contact - tunneling between the metal leads and graphene 
zig-zag edges. It can be thought of as a limiting case of strong 
coupling (a) with the tunable hybridization between contact 
and the edge. 



proximation the general problem of transmission through 
the NG interfaces and NGN structures for arbitrary over- 
lap integrals t s , t g , and t' , in square-lattice normal metal, 
graphene, and at the interface, respectively. From the 
lattice Schrodinger equation, we derive the wave-function 
matching conditions at the interfaces. We apply them 
first to the scattering off the NG interface and determine 
the reflection probability back to the normal metal. 

Then, we turn to the conductance of an NGN junction. 
It is the sum of two contributions originating from prop- 
agating and evanescent modes in graphene. The number 
of the available propagating modes is proportional to the 
radius of the Dirac cone qc taken at the gate voltage en- 
ergy eVa, 9g °c Vg- Their contribution to conductivity 
is therefore proportional to Vg', it is independent of the 
graphene length L in the limiting cases qcL <C 1 and 
qcL » 1. At the Dirac point (Vg = 0), there are no 
propagating modes; however, there is a contribution to 
conductivity from the evanescent modes that obeys the 
Ohm's law, i.e. is inversely proportional to the length 
of the graphene strip L 11 ^ 13 . For non-zero gate voltage 
Vg the contribution of evanescent modes crosses over to 
L~ 1 (qcL)~ 3 behavior for qcL ^> 1. We also observe 
that for the symmetric case t' 2 = t s t g , a narrow graphene 
strip is ideally transmitting, both through the propagat- 
ing and evanescent states. Moreover, the conductance 
of such strip precisely at the Dirac point is intependent 
of the interface overlap integral t' . Even for opaque in- 
terfaces, t' 2 <C t s t g , we find that some modes transmit 
ideally, in a direct analogy to the Fabry-Perot resonances 
in a double-barrier structure. 



II. NG INTERFACE 

We consider a zigzag interface between normal metal 
with the square lattice and graphene with the bond 

re- 



lengths a s and a g , matched such that a g 



,/V3, 



is taken near the half-filling, so that the Fermi surface is 
nearly a square, whereas the graphene is tuned close to 
the Dirac point. In the normal metal, the wave functions 
are plane waves, c(r) = (e lksX + r e~ lksX ) e lkyV , x < , 
defined on the sites of the square lattice. First and second 
terms describe incident and reflected waves, respectively. 
In graphene, one has to define two transmitted waves 
corresponding to two sublattices a and b, [d a (r),d b (r)]. 

For the zigzag interface, the wave vector component 
along the interface k y is conserved, and thus the two- 
dimensional scattering problem reduces to a collection of 
one-dimensional problems for different values of k y . The 
transverse component of wavevector (k x ) is not conserved 
and must be found from the energy conservation. 
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spectively, Fig.O We assume that the normal metal band 



FIG. 2: (a) A zigzag interface between a square and graphene 
lattices. Sites belonging to the two sublattices in graphene are 
shown by letters a and b. On the top, boundary conditions are 
sketched: We introduced two columns of fictitious amplitudes, 
ci in graphene, and d b _ 1 in the square lattice, used for wave 
functions matching at the interface, (b, c) A NGN contact of 
the width L — a g (3N + l)/2 with an odd-integer (b) or an 
even-integer (c) value of N. 



In the tight-binding model, the amplitudes obey the 
lattice Schrodinger equation. In particular, at the inter- 
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face, we obtain the following set of equations, 

Ecq = —2t s co cos k y a s — t s c^\ — t'd^ ; 



Ed a 



Ed\ /2 



-t'co — 2t g cos 



V a 9 ib 



di/2 ; 



t rl a Of rn- ^ k V a 9 ja 

" r 9«3/2 ~~ ZT g cos 2 ' 



(1) 



Near half-filling, k s a s 7r/3; however, these results ap- 
ply to the case of arbitrary filling, where k s a s ^ 7r/3. As 
anticipated, for an opaque interface t' 2 <C i a is, the re- 
flection approaches 1. Note also that one can only match 
the wave vectors in the y-direction for a very limited set 
of wave vectors around k s a s = 7r/3; all other states are 
ideally reflected by the interface. 



where we moved to discrete notations, explained in Fig. [2] 
The subscript shows the distance from the interface, mea- 
sured in corresponding lattice constants. 

Away from the interface, the solution of the 
Schrodinger equation is a superposition of plane (or 
evanescent, see below) waves, with the energies 

E s = — 2t s (cos k s a s + cos k y a s ) , 



E g - ±t g 



\ 



^/Skydq I \ZHkydq 3k„a„ 
1 + 4 cos -r— - cos -f— - + cos — ^ 



(the Dirac K-points), k g = q x , k y — ±47r/(3v / 3a 9 ) 



E± = ±3t g q/(2a g ), with q = (g~ + q y 
creteness, we only consider the "upper 

V* ^ A rrr I 



In this Article, we are mainly interested in the low-energy 
regime, E -C t g . In this regime, the energy can be ex- 
panded in the vicinity of the points where it turns to zero 

%, 

2 ) 1 I 2 . For con- 
" Dirac point, 

47r/(3v / 3a ff ); by symmetry, the "lower" one con- 
tributes to transport identically. For propagating waves, 
from the graphene (lattice) Schrodinger equation we ob- 
tain d b = ie i(p d a , where t&iup = q y /q x . 

The interface equations, Eqs. fl} can be equivalcntly 
recast in the form of the wave function matching, which 
is convenient when considering more complex scattering, 
such as NGN (see below). This can be done by continu- 
ing the wave functions across the interface. Introducing 
fictitious wave function elements c\ and d b _ 1 , the first two 
equations of Eqs. ([TJ) become 



t'd a 



t s Ci 



t'co 



These are general boundary conditions and can be ap- 
plied to a scattering problem with arbitrary arrange- 
ments. For example, scattering from N to G can be ob- 
tained if we use 
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Ak s a s 



re 



-ik s a s 



die 



ik Q a, a 



III. NGN CONTACT: PROPAGATING AND 
EVANESCENT MODES 

Consider now a graphene sheet of length L connected 
to two square-lattice metal electrodes (with the same 
overlap integral t s ) at two ideal zigzag interfaces. Note 
^2) that such an arrangement is only possible provided L = 
— a g (3N + l)/2, with an integer N. For an odd-integer 
\ N the square lattice leads are aligned (Fig. [2b) , for an 
I ■ even-integer N they are shifted by half a period (Fig. 
' [2t)- Note that sites at the two interfaces always belong 
to different graphene sublattices. We describe both cases 
on equal footing. 

We take the wave function in the left electrode in the 
same form as before, c{r) — (e lksX +re~ lk,x ) e lkyV , the 
wave function in graphene as a combination of left- and 
right-moving waves for each sublattice, and the wave 
function in the right electrode as the transmitted wave 
f(r) = w e~xjp(ik s (x — L) + ik y y). 

The equations for the NGN structure read 



t'c 
t s c\ 
t'fo 

tsf-1 



= tg{d\- 

= t'(d? - 

= tg(d? 

= t'{d\- 



-<)-! ; 
■ d")o ; 

- ^r)(3iV+l)/2 
' dr)(3iV+3)/2- 



(6) 
(7) 
(8) 
(9) 



In addition to c\ and d b _ 1 , we introduced two more "un- 
physical" amplitudes, d® 3N+3 y 2 and f-\. We use now 

d b r = —ie lip dr and d\ = ie~ %,p df, where <p in both ex- 
pressions is defined for the right-moving electrons. Solv- 
ing this scattering problem for qa g <C 1, we find for the 
transmission amplitude 



w 
J 



J- 



—2i cos (p sin k s a s e lk 

p-l e -ik.a. cos (^_ 



(10) 



(3 e s<ls cos(ip — qL cos ip) — 2 sin(qL cos ip) 



Substituting into the above boundary conditions, we find 
the equations connecting the amplitudes of the waves, 



t s {e ik ° a ° 
f(l 



■ re 
r) 



-ik s a & 



)-- 

tgd e 



= t'dg ; 



(3) 
(4) 



Thus, for the reflection coefficient R — |r| 2 , at low ener- 
gies {kgdg <C 1), we find 



Note that for qL <C 1 and t' 2 — t s t g , the junction is 
ideally transmitting. 

For a finite length of the graphene strip, there are also 
solutions which are not propagating, but rather exponen- 
tially increasing or decreasing as e ±VxX - the evanescent 
states. Their energy reads 



* 2 



R = 



+ - 2sin(^ + fe 5 a 5 ) 



P + 0- 



2 siri(ip 



Ut 



(5) 



-s ' g 



which is defined as far as |ry x | < 1^]. Similar to the 
propagating case, we can find the relation between the 
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components of the graphene wavefunction, 



±4 



1% 



q y - r\x 



(11) 



where ± correspond to E±. In the following, we con- 
sider the positive energy branch, E+. While there is no 
propagation in the evanescent case, we can still define 
the right (rj x > 0) and left (i] x < 0) components of the 
wave function - the wave "propagates" in the direction 
of decay. For these components, d b r = Zd ^, d\ = Z^ 1 df, 
where Z = ^J{q y + \Vx\)/(<ly ~ \Vx\)- The transmission 
amplitude through the evanescent modes becomes 



-2isinh£ sinfc s a s e lks<ls J 



(12) 



J = p~ 



sinh(jj x L + C) — 2 smhr] x L 



+ (3e lk ° a ° smh(r, x L-0 , 



with ( = \nZ. Clearly, for t g t s — t 12 and rj x L < 1 wc 
obtain again the perfect transmission. 



IV. CURRENT AND CONDUCTANCE 

The current through the junction is expressed as 



Ir 



eW 



(2n) 2 

e W reVa+eV 



v x \w(k s ,k y )\' 



(2n) 2 h 



eV G 



dE / dky\w{k s ,ky)\ 2 , (13) 



where W is the width of the graphene strip in the y- 
direction, and v x = Ti~ l dE s /dk s is the group veloc- 
ity. This expression includes the contribution from both 
propagating (real k s ) and evanescent (purely imaginary 
k s ) states. 

Let us first consider the contribution of propagating 
states. The integration is carried over the momenta for 
which v x > 0. In the linear regime, from Eq. (| 1 3[) we 
obtain the conductance 



G tr _ Wq G 

Gq 7T 



7T/2 



\w\ 2 cos ipdip , 



(14) 



-it/2 



with q G = 2e\V G \/ (3t g a g ) and Gq — e 2 /2itTl being the 
conductance quantum. We were able to analyze Eq. IT4]) 
analytically in the two limiting cases, short {qGL <C 1) 
and long (qcL 3> 1) junctions. 

For short junctions, q G L <C 1, the transmission coeffi- 
cient, 

2 4(3 sin 2 k s a s 



(3 1 + (3 — 2 cos k s a s ' 
does not depend on the angle tp, and we obtain 

G tr _ 2Wq G 4(3sm 2 k s a s 
Gq 7r /3 _1 + (3 — 2 cos k s a s 



(15) 



The quantity 2Wq G /n can be interpreted as a "number of 
transport channels." In the case (3=1, the conductance 
equals 2Wq G GQ/ir. 

For long junctions, q G L ^> 1, we use the fact that 
cos(q G L cos if) is a rapidly oscillating function of the an- 
gle <p. In particular, for (3 = 1, we have 



ly- 



cos 2 (psm 2 k s a s [sin 2 k s a s cos 2 <p cos 2 (q G L cos ip) 



( 1 — cos k s a a sin <p) 2 sin 2 (q G L cos ip) 



(16) 



In this situation, the integral in Eq. (|T3|) can be dis- 
cretized. Indeed, between the points (p n and <p n +i, such 
that cos</?„ = im/(q G L), the integral in Eq. (|13p can be 
easily calculated assuming that the slow functions cos ip 
and sin ip are constant and equal to cos <p n and sin <p n ev- 
erywhere except for in combination cos(q G L cos ip). Then 
Eq. (|T3|) becomes a discrete sum over the periods n of 
the function cos(q G L cos ip). Converting the sum into an 
integral (the integrand is a smooth function of n), one 
obtains 



(17) 



G tr _ Wq G (l - smk s a s ) 
Gn 4sinfc s ac cos 2 k s a s 



This result is length independent, similar to Eq. (|T5j) . 
and for k s a s = 7r/3 the conductance of a long graphene 
layer is suppressed as compared to a short layer. 

For untransparent interfaces, (3 <C 1, the transmission 
coefficient can be approximated as 



w 



A(3 2 sin k s a s cos 2 ip 



cos 2 (q G L cos ip — ip) 2 + A(3 2 sin k s a s sin (q G L cos ip) 



( 18 ) 

This expression has the structure similar to that of reso- 
nant tunneling for a double barrier: Typically the numer- 
ator is of order of (3 2 <C 1, whereas the denominator is of 
order 1, and the transmission probability is small. How- 
ever, for certain directions ip m of the wave vector, when 
the cosine in the denominator vanishes, the transmission 
becomes ideal. One can expand the expression around 
the resonance ip m to obtain the Breit-Wigner structure 
of the resonance, ip = ip m + Sip, dip <C 1, 



w 



A(3 2 sin k s a s cos 2 ip r . 



(q G Lshnp m — l) 2 Sp 2 + A(3 2 sin k s a s cos 2 tp n 



(19) 

Typically one can omit 1 as compared to q G Lsuitp m in 
the denominator. The main contribution to the current 
comes from the directions around the resonances (tp close 
to <p m ). Integrating the Breit-Wigner expressions (fT9|) 
and transforming the resulting sum over m into an inte- 
gral, we obtain 



G tr /G Q = (3Wq G sin k s a s 



(20) 



Note, that here the conductance of a long layer is pro- 
portional to (3 and thus parametrically exceeds the con- 
ductance of a short layer (proportional to (3 2 , Eq. (Tl5l) ). 
This effect is due to the resonant structure of Eq. (HI 
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Let us turn now to the contribution of the evanescent 
modes. In the vicinity of the Dirac point, the number 
of propagating states vanishes proportionally to E, and 
thus the contribution of the large number of evanescent 
states with ps \q y \ becomes dominant. At zero energy, 
E ps +0, and q y > 0, one has Z — > oo. Thus Eq. JTSJ) 
becomes 



w+ 



-2i sin k s a s e lksa " 



f3~ 1 e 



fie 1 ' 



(21) 



Interestingly, for negative values of the deviation q y < 
from the K Dirac point, the transmission amplitude for 
zero-energy states (E pa +0 and Z — > 0), is different: 



—2i sin feed, e 1 



(22) 



Although for (3 = 1 the probabilities are the same, 
|w_| 2 = |w+| 2 , for j3 7^ 1, they are different. This effect 
is related to the chirality of the surface states in a given 
Dirac valley. Indeed, a zigzag interface of graphene sup- 
ports a continuum of chiral edge states at zero energy^: 
They can only propagate in the direction q y > for the 
lower Dirac point K and in the direction q y > for the 
upper Dirac point. These states decay exponentially in 
the bulk of graphene. In a graphene nanoribbon, which 
is essentially our graphene sheet disconnected from the 
electrodes, these states become evanescent modest. It 
is this chiral nature of the surface states that cases an 
asymmetry between the two directions of q y . 

Note that the expression for w- can be obtained from 
w+ by flipping the sign of n x ; thus integration over posi- 
tive and negative values of q y in the conductivity is equiv- 
alent to integration over positive and negative values of 
rj x - The tunneling probability is 



4 sin 2 La, 



w 



p-2 e 2r, x L + p2 e -2, lx L — 2 COS 2k s d s 



(23) 



It becomes ideal for (3 — exp(r] x L). 

The conductance including now both Dirac points is 



+00 



G ev W I 2 

~Gq = 7rL ' ""' H 



W(tt - 2k s a s ) 



tan k s a s 



(24) 
(25) 



Surprisingly we find that the resulting conductivity is 
independent of the value of (3, even for small (3 (weak 
contact between the metal and graphene)! By inspec- 
tion one can see that in this limit the major contribution 
comes from W-, i.e. the negative values of q y . In the 
regime (3 <C 1, the graphene region reduces to the zig- 
zag nanoribbon, which is knowni^ to have a continuum 
of surface states for q y < with energies that scale as 
e -\q v \L _ Precisely at the Dirac point these states only 
have non-zero amplitudes at one of the sublattices. The 
conductivity is dominated by the tunneling through these 



states. A similar effect was discussed in Ref. for tun- 
neling between a carbon nanotube and a metallic elec- 
trode. 

Still, Eq. clearly cannot hold for (3 = 0. To estab- 
lish the limits of its applicability we analyze the terms 
that were dropped while going from Eq. (fT2|) to Eqs. (j2Tj) 
and (f2"2"| . Taking the limit ( ^ ±00 is only valid for ex- 
act zero-energy states. However, for finite transport or 
gate voltage such that qa <C l/L, Z(q y < 0) ps Qg/Wv\- 
The conductivity Eq. is dominated by terms with 
I Q y I ~ |??!b| ~ 1/L. The subdominant terms in Eq. (Ti"2"|) 
therefore become non-negligible when qc > (3/L. Thus, 
the /3-independent expression ([24]) only holds for low 
enough transport and gate voltages, such that qa < (3/L. 
In this regime our result agrees with the one obtained by 
Schomerus for (3 = 1—. 

To analyze the intermediate regime (3 <C qcL <C 1, 
we keep for q y < the leading term (3~ l exp(— ik s a s + 
r] x L — C) in the denominator of Eq. p2[) . This term 
only dominates provided n x is not too small, such that 
exp^L) 3> 2r) x /qG- The solution of this transcen- 
dental equation, which in the leading order becomes 
i] x = —L^ 1 ln(qaL/2), provides the cut-off in the integral 
over rj x in the expression for the conductance. Explicitly, 
we have 



16W(3 2 



Gq KL(q G L) 



■ In 2 (qcL) sin 2 k s a s 



Thus, away from zero energy, the the contribution of 
evanescent modes rapidly vanishes, and conductance be- 
comes proportional to 1 . 

For yet greater gate voltages, qcL ^> 1, the conduc- 
tance is always dominated by the states with r\ x L < 1. 
From q y — r/ x = we therefore find that at any rate 
\q y \ 3> \rj x \, and thus C = Vx/qc ^ Tj x L. This means 
that we can disregard £ in the denominator of Eq. (fl"2j) . 
and in the numerator we replace sinh £ with r\ x jqG. Writ- 
ing dk y pa ri x dr) x /qG, we find that the integrand in the 
expression for the conductance contains the third power 
of r] x multiplied with e~ n * L . Consequently, we obtain 



12((3)W 



Gq ~ TrL(q G L) 3 I/?" 1 e - lk ° a ° + (3 e lk ° a ° - 2 



(26) 



Here, £(n) is the Zeta function, with £(3) ps 1.2021. 

The results for the conductance in different regimes are 
shown in the Table. 



V. CONCLUSIONS 

We constructed the wave-function matching conditions 
at the zigzag interfaces between square (N) and graphene 
(G) lattices, and determined transport properties of the 
NGN structure, concentrating on the regimes of "ideal" 
interface (t 12 = t a t s ) and highly resistive interface (t 12 <C 
t a t s ). In accordance with earlier predictions, at the Dirac 
point the conductance is dominated by the evanescent 
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Propagating Evanescent 
= 1 qo < 1/L q G L- 1 

qg > 1/L qg L-^ggL)' 3 

(3«l?c«/)/i qaP 2 L- 1 

9G<1/L g G /3 2 f3 2 L-\q G L)~ 2 In 2 (q G L) 
go > 1/L gc/3 ^L-'igaL)- 3 

TABLE I: Functional dependence on the parameters g G , 
L, and (3 = t' 2 /(t s t g ) of contributions to the conductance 
from propagating and evanescent modes in different transport 
regimes. 

modes and scales inversely proportionally with the length 
of the contact L. However, the situation changes qual- 
itatively as soon as one departs from the Dirac point - 
for instance, by changing the electron concentration via 
the gate voltage Vg- The propagating modes start to 
contribute to the conductance. For small Vg, such that 
qoL <C 1, their contribution is length- independent; in 
particular we find that for ideal interfaces the transmis- 
sion equals one independently of the angle of incidence. 
Further yet from the Dirac point, qcL 3> 1, the conduc- 
tance is determined by the propagating modes, whereas 
the evanescent modes' contribution decays as L~ 4 . In 
particular, for qcL 3> 1 and t' 2 <C t s t g we found a regime 
similar to resonant tunneling in double barrier structures. 
As a consequence, the propagating modes' contribution 
of a long junction in this case is greater than the one for 
a short junction. 

We found that the contribution of evanescent modes 
at the Dirac point does not depend on /3, but diminishes 
as we apply the gate voltage. In the regime qcL 3> 1 this 
contribution is suppressed. 

Besides looking at the low transport bias regime con- 
sidered above, it also may be interesting to study the 



non-linear I-Vt characteristics. Our results suggest that 
with increasing transport voltage Vt, the contribution of 
the evanescent modes to the current saturates with volt- 
age beyond Vt > t g a g /L, while the propagating modes 
contribution increases as V t 2 . However, the quantitative 
discussion of non-linear regime is problematic, since the 
result would essentially depend on the potential distribu- 
tion in contacts and over the graphene sheet. Investiga- 
tion of this potential distribution would require solution 
of Poisson equation coupled to the equation for the par- 
ticle density, and goes beyond the scope of this Article. 

The zigzag interface considered in this Article is the 
simplest case of a contact: The periods of the lattices 
match, a s — V3a g , and the momentum component k y 
along the interface is conserved. In real experimental sit- 
uations both of these conditions will be difficult to real- 
ize: The interfaces are disordered, and the lattice periods 
may be incommensurate. This Article illustrates an im- 
portance of the interface contribution to the transport 
and provides the basis for future research in this direc- 
tion. 

We thank S. Trugman and H. Schomerus for useful 
discussions. The authors acknowledge Aspen Center for 
Physics, where this research was initiated. This work was 
supported in part by US DOE. 

Note added. At first glance, the result (|26[) disagrees 
with Ref. [TH, which finds G ev oc L~ x for all gate volt- 
ages. However, the chemical potentials in Ref. [l3|were ar- 
ranged in a different way than we have done it above: The 
chamical potential of graphene sheet is fixed to the Dirac 
point, whereas chamical potentials in the normal metal 
are varied. Recently, a study of non-linear transport 
by Robinson and Schomerus has been made available^. 
Both our above results and the results of Ref. [13 follow 
in appropriate limiting cases. 
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